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ABSTRACT 

By combining the recently derived X-ray luminosity function for Galactic X-ray binaries (XRBs) by Grimm 
et al. (2002) and the radio-X-ray-mass relation of accreting black holes found by Merloni et al. (2003), we 
derive predictions for the radio luminosity function and radio flux distribution (logN/logS) for XRBs. Based 
on the interpretation that the radio-X-ray-mass relation is an expression of an underlying relation between jet 
power and nuclear radio luminosity, we derive the kinetic luminosity function for Galactic black hole jets, up 
to a normalization constant in jet power. We present estimates for this constant on the basis of known ratios 
of jet power to core flux for AGN jets and available limits for individual XRBs. We find that, if XRB jets do 
indeed fall on the same radio flux-kinetic power relation as AGN jets, the estimated mean kinetic luminosity of 
typical low/hard state jets is of the order of ( Wxrb )~ 2 x 10 37 ergs with a total integrated power output 
of Wxrb ~ 5.5 x 10 38 ergss~ 1 . We find that the power carried in transient jets should be of comparable 
magnitude to that carried in low/hard state jets. Including neutron star systems increases this estimate to 
HxRB,tot ~ 9 x 10 38 ergs s _1 . We estimate the total kinetic energy output from low/hard state jets over the 
history of the Galaxy to be £xrb ~7x 10 56 ergs. 
Subject headings: 



l. introduction 

X-ray binaries (XRBs) have long been thought of as clas- 
sic examples of accretion flows with little or no complicated 
outflow physics in the way of understanding the dynamics of 
these systems. It was not until recently that jets and outflows 
were realized as important ingredients in the process of accre- 
tion even in XRBs. The reason for this late appreciation is the 
fact that black hole XRBs ar e less radio loud than A GNs by 
about 4 orders of magnitude (Heinz & Sunvaev 2003). 

Accreting neutron star binaries, whic h are now also known 
to regularly produce radio emitting jets ( Fen derl2005l) . are yet 
anothe r factor of 30 more radio quiet than bla ck hole XRBs 
jFender & Kuulkersll200!l iMigliari et al]l200l . This further 
reduction in radio flux has made quantitative analysis of neu- 
tron star jets difficult to the point that rather little is known 
about them at present. Therefore, we will henceforth focus on 
the properties of jets from black hole XRBs throughout most 
of this paper and will present an extension of our work in i|3.5l 
to neutron star XRBs. 

In consequence, while XRBs are among the brightest X-ray 
sources in the sky, their radio output is rather weak. This, and 
the fact that the number of Galactic XRBs is much smaller 
than the number of AGNs, are to blame for the fact that very 
little information about the radio power from XRB jets exists 
and that the physical properties of these jets are even less well 
known that those of AGN jets (which by themselves still pose 
many questions as to their exact makeup and dynamics). 

Two types of radio loud jet-like outflows have been ob- 
served in black hole XRBs. Following transitions in the 
X-ray spectral state, one can often observe bright radio 
flares, which are optically thin and show a powerlaw tem- 
poral decay on day timescales. The radio emission can 
be re solved on sub-arcsecond scales and shows proper mo- 
tion (Mirabel & Rodriguez! 1 19941 iHiellming & Rupenlfi995t 
iFender et alll999l) . These ejections can carry a large amount 



of power, but are relatively rare. The second kind of jet- 
like outflow is associate d with the so-called low/hard s tate 
in black hole XRBs (see McClinto ck"& Remillardl2003l for 
a thorough review on X-ray states in XRBs). These outflows 
show a flat, optically thick radio spectrum and typically do 
not show a temporal decay. Only two of these sources have 
been r esolved, showing a collimated jet on mi l li-arcsecond 
scales Stirling etalJl2001l iDhawan et ai]l2000t iFuchs et alJ 
2003). Consequently, these types of flows are referred to as 
"steady/compact" jets, and for the remainder of this paper, we 
will be mostly concerned with these types of jets. 

Steady, compact, flat spectrum jets are observed in virtually 
every Galactic black hole XRBs in the low/hard state that is 
accessible to radio instruments. Significant progress in our 
understanding of XRB jets has recently been made when a 
relation between the radio luminosity from these steady, com- 
pact jets and the har d X-ray luminosity was discovered in 
low/hard state source dCorbel et alJl2002t iGallo et alJl2003l) . 
This relation has since been found to extend from stellar mass 
blac k holes all the way up the b lack hole mass scale to AGN 
jets tMerloni et al. 2003; Falcke et al 120041) . In the remainder 
of this paper we will refer to this relation as the "fundamental 
plane of black hole activity" (FP hereafter). 

The FP relation expresses both the fact that the radio lu- 
minosity from compact jets is non-linearly correlated with 
the X-ray flux, as had already been discovered in the case 
of XRBs, as well as the non-linear dependence of radio flux 
on black hole mass. The latter is responsible for the fact that 
AGN are so much more radio loud than XRBs. These rela- 
tions can be understood naturally if the physics underlying 
jet formation (ultimately, strong gravity and MHD in the in- 
ner ac cretion flow) are invaria nt under changes in black hole 
mass (Heinz & Sunvaev 2003). 

This scale invariance implies that a relation between the ra- 
dio luminosity L r emitted by the jet and the kinetic jet power 
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Wkm is underlying the radio-X-ray-mass relation. For rea- 
sonable parameters 1 , this relation takes on the form L ra d oc 

Wh^~ a ' , where a T ~ is the radio spectral index. To- 
gether with the recent discovery of the FP, this relation can 
provide a powerful diagnostic for studying the properties of 
jets. 

A second recent insight that has spurred heightened interest 
in the statistical study of XRBs is the firm det ermination of 
the M ilky Way XRB X-ray luminosity function ( Grim m et al.l 
120021) . In this paper, we will couple the information contained 
in the XRB luminosity function with the predictive power of 
the FP relation to derive some of the missing statistical proper- 
ties of the XRB jet population needed in planning observation 
campaigns for XRBs and in modeling the jets of XRBs and 
their impact on the interstellar medium. 

In ij2]we will derive the kinetic luminosity function for 
Milky Way low/hard state XRB jets, parameterized by the un- 
known normalization of the kinetic power - radio power rela- 
tion underlying the fundamental plane relation. In S|3] we will 
present an estimate of that normalization derived from AGN 
jets and the available limits for individual XRBs, and apply it 
to the XRB jet population to derive the absolute kinetic lumi- 
nosity function. In 2]we derive, as corollaries to the kinetic 
luminosity function, the predicted radio luminosity function 
of XRB jets and the predicted Galactic radio logN/logS dis- 
tribution. Section|5]summarizes our results. 



2. DERIVING THE KINETIC LUMINOSITY FUNCTION OF BLACK 
HOLE XRBS IN THE LOW/HARD STATE 

The discovery of the FP relation has enabled a num- 
ber of diag nostic tools to be d eveloped for accreting 
black holes (lMerlonill2004t lMaccarondl2004t iHeinzl 120041 
iHeinz & Merlonil2004HMarkoffl2005l) . Since the current the- 
oretical understanding of the FP relation is based on a rela- 
tion between kinetic jet power and radio luminosity, it is now 
possible to convert a radio luminosity f unction into a kinetic 
luminosity function (Heinz et al. 2005). This is a powerful 
and useful technique for AGN jets, where we can measure 
the radio luminosity function and have good estimates of the 
kinetic power in a number of sources. In XRBs we have nei- 
ther, due to the characteristically lower radio fluxes of XRB 
jets (which has delayed their mass discovery by 30 years rel- 
ative to the large known sample of AGN jets) and the lack of 
good estimators of kinetic power in these sources. 

However, the small dispersion in black hole mass in black 
hole XRBs and t he recentl y published X-ray luminosity func- 
tion for XRBs (Grimm et al. 2002) allow us to use the FP re- 
lation to construct both a predicted radio luminosity function 
and a kinetic luminosity function. These estimates are nec- 
essarily subject to considerable uncertainty from a number of 
unknowns: (a) The conversion of the total XRB luminosity 
function to a black hole luminosity function (essentially, the 
black hole fraction of XRBs), (b) the normalization of the ki- 
netic jet power for a given radio power (this holds for AGNs as 
well), and (c) the mass distribution of stellar mass black holes. 
We will parameterize our results in terms of the fiducial values 
we assume for these unknowns to allow for easy adjustment 
with future improvements in the knowledge of these parame- 
ters. 

1 For an electron spectrum with powerlaw index p = 2 and proportionality 
between kinetic jet power and magnetic energy density, W^in oc B 2 , i.e., for 
magnetically driven jets 



2.1. The XRB luminosity function 
The FP relation implies that below a certain critical accre- 
tion rate m cr it = M cr it/MEdd ~ 0.01, there exists a correla- 
tion between the mass of a black hole M, its radio luminosity 
L T , and its X-ray luminosity L K of the form 
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L T = L l^M 
L = 1.Q x lO^ergss" 1 



(1) 

(2) 



where, for convenience, we expressed the X-ray lumi- 
nosity in terms of the Eddington luminosity LEdd = 
1.3 x 10 38 ergss _1 of a one solar mass o bject: lx = 
LyilL v M. Following the interpretation of iMerloni et alJ 
(2003); Fe nder et alJ ( l2004bl) . we will assume in the follow- 
ing that the transition at m cr i t is due to a state change as ob- 
served in XRBs, where low luminosity sources switch into 
the low/hard state. In XRBs, which we are concerned with 
here, steady, flat spectrum jets are closely associated with the 
low/hard state. The details of this transition are irrelevant, 
we will simply use the associated X-ray luminosity as an up- 
per bound on the luminosity function. For a fixed black hole 
mass M — 10 Mq M\q, this transition corre sponds to a criti- 
cal X- r ay luminosity l-x r .r\t = 10 Mi n m„\+. iMivamot o et al.l 
( 119951) : iMaccarona ( 120031) showed that the transition from 
high/soft to low/hard state actually occurs over a range of val- 
ues for TO cr it, and that a hysteresis exists with sources on the 
descending luminosity branch staying in the high/soft state 
longer (to lower luminosities) than ascending sources, which 
mak e the state transition fr om hard to soft at higher luminosi- 
ties. iFender et alJ ll2004bl) show how the properties of XRB 
jets correlate with the position of a source on this hysteresis 
track. We will comm ent o n the impact the uncertainty of m cr j t 
has on our re sults in £13.31 

Following IHeinz & Sunvaev) (120031) . the radio luminosity 
L T of the flat spectrum core of the jet is related to the jet power 
^jetby 



jet 
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where a T = — <9(ln_L„)/9(ln^) ~ is the radio spectral in- 
dex an close to zero for the flat spectrum sources under con- 
sideration. Henceforth we will use a r = unless noted other- 
wise. Wq is currently not well known: it is the normalization 
of the kinetic jet power in relation to its radio luminosity 2 . 
It is not well known because jet power is still a very diffi- 
cult quantity to measure, even after four decades of research 
on jets. We will attempt to estimate it below and will carry 
it through the algebra as a parameter until then. Combining 
eqs. Q and y}, we can now write 
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jet 
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M 0.6 



(4) 



We will use the X-ray luminosity f unctio n for the Milky 
Way XRBs provided by Gri mm et alJ (|2002) (see Fig. 12 in 
their original paper). Where appropriate, we provide esti- 
mates based both on the actual data they derived (numerically 

2 Naively, Wo can be taken as a radiative efficiency, though it should be 
kept in mind that most of the radiation is emitted at high frequencies, thus the 
radiative efficiency is often times not well defined, as the high energy cutoff 
of individual sources is hard to measure and it is not clear which parts of the 
spectrum actually arise in the jet and which arise in the accretion flow. 



The kinetic power of XRB jets 
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integrated) and parameterized using the powerlaw approxima- 
tion dN/dl = NqI' 13 , fitted bv lGrimm etafl <2002li to give: 



dN HMXB 

dl x 

dN^MXB 

dL 



1.6 



(5) 
(6) 



for l x < 1 as well as the actual data. The X-ray luminos- 
ity function does not show a low luminosity cutoff. But at 
low luminosities the observed number of sources is very small 
due to the sensitivity limits of the ASM on RXTE. However, 
a comparison with higher sensitivity (but sma ller sky cover- 
age) A SCA data of the Galactic Ridge Survey ( Sugiza kTet al.l 
120011) show that the XR LF does not req uire a cutoff down to 
~few 10 33 ergs s" 1 (iGrimm et all2002l) . 

The radio-X-ray correlation in XRBs has been observed 
down to X-ray luminosities of about l x > 10—4 in the case 
of GX 339-4. While there is reason to believe that the radio- 
X-ray relation might break down at very low luminosities (the 
synchrotron X-ray component from the jet scales more slowly 
with m, F x sync h oc m 1-18 than the X-rays from the accre- 
tion flow, Fx.acc o c m 2 ~ 2 3 lMerloni et al.ll2003l lHeinzll2004l 
I Yuan & Cul 120051) . we will take l x ,min = 10~ 4 as a secure 
upper limit on the low luminosity cutoff of the radio-X-ray 
relation. 

The relative fraction of black holes as a function of l x 
in both the LMXB and HMXB luminosity functions is un- 
known. Thus, for lack of better knowledge, we will use the 
ratio of observed black hole to neutron star XRBs, which 
is about 10%. We parameterize this unknown fraction as 
( = [cOVbh / dl x ] I '[cL/Vns / dlx] = 0.1£o.i and assume that the 
fraction of black holes is constant as a function of l K , Given 
the uncertainty in these estimates and the lack of knowledge 
of the black hole mass function, we will henceforth assume 
that the mean black hole mass is 10 MqMiq and evaluate all 
quantities for this black hole mass. 

The instantaneous kinetic luminosity function for jets from 
Galactic black holes is then 
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for l^WoM 55 < Wj et < m° r f W M°- 55 . For the power- 
law parameterizations, we have 
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in the same bounds. All that is left to do is the determination 
of Wo — clearly the most difficult and uncertain part of this 
exercise. 

3. ESTIMATING THE ABSOLUTE KINETIC LUMINOSITY 
FUNCTION 

3.1. Estimating Wo 
3.1.1. Estimates of the kinetic power in AGN jets 

Information on the kinetic power from XRB jets is just 
starting to become available. We will discuss the av ailable 
estimates and how they can be used to constrain Wo in A3. 1.21 
First, however, we will consider the case of AGN jets, for 



which information about the kinetic power is readily avail- 
able and undisputed in a number of important cases. There 
are three AGN j et sources in th e sample used to derive the FP 
relation by Merloni et al. ( 2003) which have reliable estimates 
of the jet power: M87, Cygnus A, and Perseus A. These esti- 
mates are derived from average powers on very large spatial 
scales in the cases of Cygnus A and Perseus A, while for M87 
power estimates exist both for larg e spatial scales fr om X-ray 
cavities, and for jet scales (Bicknell & Begelman 1996). 

We can use these sources to estimate Wo ( Hein z et"al] 
120051) . We should keep in mind, though, that the sources 
each have some offset from the fundamental plane relation 
due to the intrinsic scatter around the plane. Furthermore, 
the large scale power estimates are time averaged (over times 
much longer than the dynamical time of the central jet that 
is currently producing the radio emission) and might not be 
representative of the current jet power. However, the fact that 
the estimate we get from M87 alone (where the large scale 
power estimate agrees with the smaller scale estimate) is con- 
sistent with the ones obtained from the other two sources in- 
dicates that we are probably not too far off. It should also be 
noted that all of these sources show compact, flat spectrum 
jets, which is critical if we want to use them as templates for 
low/hard state jets from Galactic XRBs. 

The rough power estimates for the three sources are: 
W ict » 10 45 - 7 ergss-\ lO^ergss" 1 , an d 10 44 - 5 ergs s" 1 
for Cy g nus A, M87, and Per A respec ti vely (ICarilli & Barthell 
1 19961 iBicknell & Begelmanl 119961 iForman et all 120051 
Fabia n TeTaiTl2002l) . while their core radio luminosities are 
10 414 , 10 39 8 , and 10 417 ergs s _1 , respectively. Thus, 
the normalization constants we derive for the sources are 
Wo.CygA ~ 6.3 x 10 37 ergs s" 1 , W QM S7 ~ 2x 10 37 ergs s" 1 , 
and Wo.PerA ~ 2.6 x 10 36 ergss _1 , each carrying consider- 
able uncertainty from the estimate of the kinetic power and 
from the fact that the sources scatter around the fundamental 
plane. 

We can also calculate what the radio power should be if 
there were no scatter around the relation, by estimating the 
unbeamed radio flux from the X-ray luminosity and the mass 
of the objects, using the FP relation. In this case, the normal- 
ization constants are Wo,cygA ~ 10 38 ergss _1 , Wq 7 ms7 ~ 
6.3 x 10 37 ergss-\ and W ,p C rA ~ 3.7 x 10 37 ergss -1 , 
which is significantly more homogeneous. 

This expression is more appropriate to use in an average 
sense when considering a large sample of sources, since the 
radio emission probably carries the biggest fraction of the 
scatter in the fundamental plane relation due to variations in 
relativistic beaming, black hole spin, and spectral index. In 
this particular case it is also more appropriate because we 
are using the X-ray luminosity of XRBs and the fundamental 
plane to convert X-ray fluxes into radio fluxes. Thus, using X- 
ray fluxes to calibrate this method seems to us to be the most 
consistent approach. We will therefore use this "corrected" 
estimate of the efficiency, which averages to 

W w 6.2 x 10 37 ergs s" 1 W 37 . 8 (10) 

with about an order of magnitude uncertainty. We will carry 
W37.8 through the rest of the paper to allow the reader to 
adjust for future improvements and different personal pref- 
erences in this value. 

3.1.2. Limits from power estimates of XRB jets 

There is relatively little information about the kinetic power 
output of individual Galactic XRBs that we could use to cal- 
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ibrate Wq. The best information available to date is on the 
jet in Cygnus X-l. The compact jet is resolved on scales 
of 15 mas, and using a si mple Blandford-Koenigl model 
jBlandford & Koenigllll97 9) to describe the radio emission, 
one can put a lower limit of Wj e t > 3 x 10 33 ergs s -1 , which 
is very likely much lower than the true kinetic power. For 
the radio flux of 12 mJy, this trans lates to a lower limit of 
Wo > lO^ergss" 1 Stirling et all2001l) . 

Based o n the IR observa t ions o f compact jet in GRS 
1915+105. iFender & Hendrvl feOOOl) argue that the normal- 
ization for the kinetic power must be larger than Wq > 2 x 
10 35 ergs s _1 , since, averaged over a sufficient time scale, the 
kinetic power must exceed the radiative output of the jet. Sim- 
ilar argum ents are commonly used in the analysis of Blazar 
emission jGhisellini & Celottl2001l) . 

Detailed modelin g of the spectral pro perties of th e XTE 
J 1 1 1 8 jet has led Markoff et1ril(l200~il) and lYuan & Cull (120051) 
to propose values of the kinetic power that translate to normal- 
ization values of Wq fss 10 38 ergs s _1 and 2 x 10 37 ergs s _1 , 
respectively. It is important to note that the jet power in 
these models is essentially a free parameter, since param- 
eters like the proton content and the opening angle of the 
jets are unkown. Nonetheless, these estimates show that 
reasonable parametrizations used in the literature for a vari- 
ety of mode ls are consistent with the estimate we presented 
in eq. iTOI . iMalzac et al.1 (2004) argue for a lower limit of 
W) ^ 2 x 10 38 ergss _1 [higher but still consistent with our 
estimate in eq. (II Oi below] on the basis of their model for the 
timing behavior of the accretion flow. However, this model 
does not differentiate between a wind-like outflow and a jet 
(as long as the flow acts as an energy sink) and it is not clear 
how this limit can be turned into a constraint on the jet power 
alone. 

Calorimetric observations (i.e., radio lobes, X-ray cavities, 
and shocks surrounding the cocoons of radio sources) are the 
most reliable way to estimate jet powers. Recently, evidence 
for the interaction of the Cygnus X-l jet with its environment 
was discovered in the form of a ring of thermal emission that 
is probably the shock driven into the interstellar m edium by a 
so far undetected radio lobe around Cygnus X-l (Galloet al. 
2005). The analysis of this shock indicates that the aver- 
age jet power is 10 36 ergs s" 1 < (W) < 10 37 ergss 



As we will show in £13 .21 this corresponds to a range of 
5 x 10 36 ergss _1 < Wo < 5 x 10 37 , which is consistent 
with estimate we present in eq. (1101 . 

While Circinus X-l does appear to possess a radio lobe, it 
does not seem to show a similar shock as Cygnus X-l does, 
and thus estimates of the source age are more difficult. Fur- 
thermore, it is still not entirely clear whether Circinus X-l is a 
neutron star or black hole. JHeinzl2002l) estimated the source 
power to be (W) > 10 35 ergss _1 . This limit is consistent 
with that on Cygnus X-l if the system is in fact a neutron star 
rather than a black hole, in which case we would expect the 
average source power to be smaller by about the mass ratio, 
roughly an order of magnitude. 

Finally, ( Kais er et alJ 120041) estimate the mean kinetic 
power of the jets in GRS 1915+105 from the tent ative associ- 
ation with what appear to be two IRAS hot spots (Chatv et al. 
1200 ll note that this association leads to distance estimates that 
are in marginally inconsistent with other measurements) to 
be 10 36 ergs s _1 (W) < 10 37 ergss _1 . If a significant po- 
tion of the power in the GRS 1915+105 jet come from the 
compact/steady flow, this estimate is consistent with the es- 
timate from above for Cygnus X-l and again implies that 
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FIG. 1. — Estimate of the cumulative kinetic luminosity function for jets 
from low/hard state black hole XRBs for HMXBs (thin black line), LMXBs 
(thin grey line) and the total (thick grey line), and powerlaw approximation 
for high luminosities. Also shown (top axis) are the corresponding radio 
luminosities at 5 GHz, for which this plot shows the cumulative radio lumi- 
nosity function. 



5 x 10 36 ergss _1 < Wo < 5 x 10 37 ergss _1 , and it is also 
very much in line with what we will estimate below based on 
AGN jets. 

3.2. The absolute kinetic power output of Galactic black 
holes in the low/hard state 

With the estimate of Wo in hand, we can evaluate a number 
of important quantities. The first is simply the kinetic lumi- 
nosity function for XRBs. The approximate analytic expres- 
sion is given by eqn. |7} with the value of Wo from eq. Ill 0b . 
Using the actual d ata of the XRB luminosity function by 
Grimm et al. (2002), we can also plot the kinetic luminosity 
function, as shown in Fig.^ 

The second is the total kinetic power released into the 
Galaxy by steady, flat spectrum jets from XRBs. It is sim- 
ply the integral over the luminosity function: 



W tot = J dW ]C 



dN 



dW; 



jot 



Whmxb : 
Wlmxb ; 



: 2 x 10 38 ergs > 
' 3.5 x 10 38 



War.sCo.iA^o 55 (11) 
ergss^War.sCo.iA/iV 5 (12) 



Note that the power output from HMXBs is dominated by 
low luminosity sources - the total power estimate depends 
on the assumed lower cutoff and is thus a lower limit. For 
LMXBs, the total power is essentially a logarithmic function 
of the lower limit on W to lowest order, i.e., each decade in 
Wkm contributes the same amount of total power. At low lu- 
minosities, however, both the HMXB and the LMXB lumi- 
nosity functions turn over, the radio-X-ray relation presum- 
ably breaks down, and the total number of XRBs is limited, 
so, as expected, the kinetic luminosity function does not di- 
verge. 

It is worth pointing out that the value for Wtot is somewhat 
larger than the estimated t otal energy output from l ow/hard 
state XRB jets presented in Heinz & Sunvaev (2002), which 
was derived from very different arguments (we should note 
that the estimate presented in this paper should be signifi- 
cantly more reliable). This also implies that the statements 
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made in that paper about the total contribution from XRB jets 
to the Galactic cosmic ray spectrum still hold (i.e., that the 
total contribution is small, however, spectrally unusual signa- 
tures like narrow lines or Maxwellian features will stand out 
measurably). 

Given the estimated star formation rate of 3M Q yr _1 of 
the Galaxy, and the fact that the HMXB luminos ity function 
is a g ood estimator of the star formation rate ( Grimm et al. 
2003), we can estimate the total kinetic power output in other 
galaxies with known star formation rates by multiplying the 
HMXB kinetic luminosity estimate from (II Q by the ratio of 
star formation rates relative to the Milky Way. The LMXB 
luminosity function does not scale like the star formation rate. 
Instead, it should scale like the total mass in stars (dominated 
by low mass stars), of which the LMXBs are a given fraction, 
which is slowly increasing with time beyond a stellar age of 
about a billion years. 

Using the stellar mass of the Milky Way, we can also cal- 
culate the time integrated total energy that was released by 
low/hard state jets: The total stellar mass in the Milky Way 
is about 10 11 MqAAu. For a current star formation rate of 
A4 w 3Af Q yr _1 we can multiply the current power esti- 
mate from HMXBs by the star formation age of the Milky 
Way to derive an estimate of the total energy released by 
HMXB jets throughout the history of the Galaxy: 

£ H MXB,tot « 2.2 x lO^ergsMu/MsWsrsCo.iM^ 55 

(13) 

while for LMXBs, we should multiply the current energy 
output by the age of the LMXB population, which is about 
10 10 tio yrs such that 

£ L MXB,tot - 2.3 x 10 56 ergsi 10 W 37 . 8 Co.iMi°o 55 (14) 
The third quantity worth estimating from this distribution is 
the mean jet power of individual XRBs. The main uncertainty 
here is the variability of individual sources. Since the XRLF is 
a snapshot of the X-ray output from a large sample of sources, 
we will have to make an assumption about how variability of 
individual sources factors into the XRLF. For lack of better 
knowledge and for reasons of simplicity, we will consider two 
possible, simplified scenarios: 

• The XRLF reflects the temporal X-ray luminosity dis- 
tribution of each individual source. I.e., each source 
spends a fraction of its life proportional to dN/dl x (l x ) 
in the luminosity bin [l x ,l x + dl x }. The luminosity 
range spanned by transients is certainly comparable if 
not larger than the range spanned by l m i n and m cr it, 
supporting this picture. Since the XRLF is rather steep, 
in this scenario an individual source therefore spends 
most of its life at low luminosities - as has long been 
known in the X-ray community to be the case for tran- 
sient sources. 

The average kinetic luminosity of an individual source 
is simply the total kinetic luminosity of all XRBs di- 
vided by the number of black hole XRBs in the lumi- 
nosity interval 10 -4 < l x < 10m cr ; t Mio: 

(WOhmxb « 1.7 x HFergss" 1 W 37 . S M° 55 (15) 
and 

(WOlmxb w2.8 x 10 37 ergss- 1 W37.8Afi6 88 (16) 

These are very large numbers indeed, about 1-2% of 
the Eddington luminosity for a 10 solar mass black 



hole, if our estimate of W37.8 is at least of the right 
order of magnitude. This indicates that the kinetic en- 
ergy transported by XRB jets is comparable to or larger 
than the radiative output: the integral over the lumi- 
nosity functions in eqs. (0 and (|6j gives average X- 
ray luminosities of (£ x ,hmxb) ~ 7 x 10 35 ergs and 
(£x,hmxb) ~ 1-2 x 10 37 ergs. It also implies that the 
kinetic energy released in low luminosity states is at 
least comparab le to the energy a dvected into the black 
hole (see also Fen der et all20 03). 

We can now compare these numbers with the lim- 
its from Cygnus X-l and GRS 1915+105 from the 
previous section: (W)cygX-i > 10 36 ergss _1 and 
10 36 ergss _1 (H / )i9i5 < few x 10 37 ergss~ 1 , which 
imply 0.1 < W37.8 < few, very much consistent with 
the estimate derived from AGN jets. 

In this scenario, all sources in either the HMXB or 
LMXB class are assumed to be very similar, which, 
given the diversity of X-ray binaries, is clearly a sim- 
plification. The fact that the XR LF is not altered by 
the variability of individual source (Grimm et al. 2005) 
somewhat supports this simplified picture and the as- 
sumption that we can average over the luminosity func- 
tion to remove the effects of variability. 

• Alternatively, we could consider a scenario where each 
source varies around an interval in Wkm much narrower 
than the width spanned by l m i n and m cv it. In this case, 
individual sources differ from each other in average ki- 
netic power, and the distribution of (W) is identical to 
the distribution of instantaneous power W . Still, the en- 
semble averaged kinetic luminosity function $((Wj et )) 
is identical to the snapshot kinetic luminosity function 
^(Wkin) an d the number derived for the total kinetic 
power is identical to eq. (II It . 

Neither of these scenarios is likely to be entirely correct; 
however, the range in Wkin spanned by the kinetic luminosity 
function in eq. (0 is relatively narrow and it is rather likely 
that individual sources traverse a range in W- ]C t that is at least 
as large as this. Thus, using the average value from eq. J 1 5i 
for individual sources to estimate the effects on the interstellar 
medium appears to us to be a reasonable choice. 

Finally we stress again that these estimates are based on 
an extrapolation from AGN jet power estimates, which them- 
selves are somewhat uncertain and are thus working estimates 
only. Obviously, a much better and more accurate way to esti- 
mate Wkin for XRBs is a direct determination. We anticipate 
that conclusiv e measurements of Wq from XRB radio lobes 
dHeinzll2002l) are going to become available in the near fu- 
ture. 

3.3. Hysteresis and the critical accretion rate 

As hinted at in §1] the exact location of m C nt is not only 
uncertain, sources actually show a hysteresis in transitioning 
from the high/soft state to the low/hard state and back. The 
transition from high/soft to low/hard typically occurs at lower 
luminosities than the reverse transition. Before elaborating as 
to how much this will affects our estimates, it is worth exam- 
ining quickly what the effects of different values for m cr it on 
(W) and Wtot actually are. 

We have evaluated eqs. (II It and J12I for a range of values 
of rhcrit and fitted the results with a quadratic in log-log space. 
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For the total jet power, we find that 

log Whmxb w 38.34 + 0.05 [log(m crit ) + 2] (17) 
-0.238 [log(m crit ) + 2] 2 + log W 37 . 8 ( .iM° 55 

log W L mxb w 38.58 + 0.38 [log(m crit ) + 2] (18) 

-0.151 [log(m crit ) + 2] 2 + log W 37 . 8 Co.iMio 55 

Since the mean jet power is simply the total power divided by 
the number of binaries, it has the same dependence on m cr it, 
but zero order different normalizations, as given by eqs. J15i 
and 02). 

This demonstrates that our estimates are not very strongly 
dependent on rh cx it for values of m cr it > 0.01. For LMXBs, 
the power changes by a factor of 1.7 when increasing m cr j t 
from 0.01 to 0.1. Given that the expected range in which this 
hystere sis occurs is about 0.01 < m cr it < 0.1 ( IFender et al.l 
2004b), we can estimate the effect by assuming that sources 
spend an equal amount of time on both the high/soft and the 
low/hard branches. Thus, a crude estimate would imply that 
about 50% of the sources above m cr jt > 0.01 are still in the 
low/hard state and produce jets. This would increase our es- 
timates of the jet power for LMXBs by 36%, well within the 
uncertainties of our estimates. HMXBs are essentially un- 
affected by this change because there are very few sources 
above m = 0.01. 

3.4. Transient Jets 

As mentioned in the introduction, some XRBs display tran- 
sient jet emission associated with rapid changes in their X-ray 
luminosity. These events can reach radio fluxes that are much 
larger than in the case of steady, compact jets produced in the 
low/hard state. T he classic example of this k i nd of source 
is GR S1915+105 (iMirabel & Rodrfgue3 H994l IFender etaD 
1999), the complex behavior of whic h is too diverse to be re- 
viewed in the scope of this paper (see Fender & Belloni 2004, 
for a review of GRS 1915+105). Other transient sources for 
which this beha vior has be en observed include GR O J 1655-4 
jHiellming & Rupenlll995l). V4641 Ser llOrosz et all 120011) . 
XTE J 1748-288 JFender & Kuulkersl2001l). and XTE J1550- 
llCorbel et al.l20 02) (see table 1 in lFenderl200 5, for more de- 
tails). The radio emission from these events is typically opti- 
cally thin and can be spatially resolved. 

In seve ral cases, superluminal proper motion has been 
detec t ed (IMirabe l & Rodriguez! 1 1 9941 (HjeJlming^Rupen 
1 19951 IFender et alJ 1 19991 ICorbel et all 120021 iFender et al.1 
2004a), which implies relativistic velocities with Lorentz fac- 
tors r ss 2 — 5. The relative abundance of measurable pa- 
rameters in the optically thin, partially resolved case makes it 
possible to estimate the total energy contained in the emission 
region and estimating the kinetic power in these transient jets. 
We shall briefly discuss these estimates and their implication 
for the total jet power from XRBs. 

It is currently not understood what the relationship between 
the transient jet events observed during state transitions, and 
the steady compact jets observed in the low/hard state is. 
It is important to stress that we cannot use the power es- 
timates from the transient events to normalize the low/hard 
state kinetic luminosity function. Because we have limited 
the kinetic luminosity function derived above to the luminos- 
ity range typically spanned by the low/hard state, our kinetic 
luminosity function and the quantities derived from it don't 
account for the separate component contributed by these tran- 
sient jets by construction. 



There are several reasons why we cannot simply expand our 
treatment to include these sources: (a) By their nature, these 
events are transient, and unlike in the low/hard state, a source 
may or may not be emitting a transient jet at a given X-ray 
luminosity in a given state. The fraction of time a transient 
source is emitting a transient jet at a given X-ray luminosity 
is unknown, (b) The radio-X-ray relation that we used to de- 
rive the relation between the X-ray luminosity and the radio 
luminosity does not hold for transient jets, (c) The emission 
from transient jets is optically thin, and thus the relation be- 
tween radio emission and kinetic power in eq. does not 
hold (which is one of the reasons why the radio-X-ray rela- 
tion breaks down). 

In the absence of statistical information about the relative 
duration, brightness, and kinetic power of individual transient 
jet events (which would be necessary to calculate the inte- 
grated transient jet power), we can still make an educated, 
though very rough guess of the contrib ution of these jets to 
the total kinetic power: For one thing, (Fender 2005) show 
that sources undergoing a hard-soft transition are known to 
produce transient jets, while they speculate that sources mak- 
ing a soft-hard transition do not. Thus, over some luminosity 
range at the transition luminosity from hard to soft state, we 
should expect a transient jet to be emitted in between 50% to 
100% of the sources . 

Secondly, Fender (2005) show that a relation between X- 
ray power and kientic power exists in transient jet sources and 
can be (roughly) calibrated to give 

W ict = 1.3ZW°- 5±0 - 2 (19) 

This relation is surprisingly close to the relation we 
used above for the low/hard state jets: Wj e t = 
1.8L F M l° A2 Wz7 sMfn 55 , as already pointed out bv lFenderl 
(2005). 

Thus, if we knew the duty cycle £ of transient events (i.e., 
the fraction of source above m cr it which are producing tran- 
sient jets), we could graft these two functions together at m cr ;t 
to make one continuous kinetic luminosity function and inte- 
grate to get the total power. Since we have neither a functional 
form of the distribution of m CI i t nor the duty cycle of transient 
events above m cr it, we can only state that the expression in 
eq. dl9l is statistically identical to eq. and that the contri- 
bution from transient jets to the kinetic power should therefore 
be well described by eqs. (I17> and Jl 81 multiplied by £, which 
describe the dependence of the kinetic power on the value of 

TO cr it. 

As already pointed out in £13.31 increasing m cr it by an or- 
der of magnitude only increases the kinetic power estimates 
by a factor of ss 1.7. Since the transient jet source fall on 
essentially the same kinetic luminosity function with an un- 
known duty cycle £ < 1, the ratio of the kinetic power carried 
in transient jets to that carried in compact, steady jets should 
be roughly 0.7£/W37.g. This gives a total power estimate of 
order 4£ x 10 38 ergs s -1 for transient jets. 

On the other hand, we know from the observed num- 
ber of outbursts that the total kinetic power from transient 
jets in the Galaxy is of the order of few x 10 38 ergss _1 
llHeinz & Sun vaevl 12002). compared to the total estimated 
power of w 5 x 10 38 ergss -1 . Thus, based on the estimates 
of Wo we provided abov e and the kinetic po wer normaliza- 
tion of transient jets by Fend er et alJ y004b), the integrated 
kinetic powers of transient jets and of compact steady jets are 
of the same order of magnitude and the duty cycle of transient 
jets should not be much smaller than S3 10%. 
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3.5. Neutron Stars 

As mentioned in ^ neutron stars are even less radio loud 
than black hole XRBs at typical X-ray luminosities. At the 
current time it is not entirely clear whether the FP relation of 
eq. Q is applicable to neutron star systems as well: it is possi- 
ble that a radio- X-ray relation exists fo r neutron stars, but has 
a different slope (Migliari et alJl2003l) . it is also possible that 
no such relation exists at all. In either case, an extension of 
our method is not easily possible. For simplicity, however, we 
will assume that the same correlation shown in eq. Q holds 
for neutron stars as well, but with a different normalization. 

Given the difference in radio loudness of a factor of ~ 30 
at a fixed X-ray luminosity and taking into account the mass 
difference between neutron stars and black holes of about 
AIbhJM~ns ~ 10/1.4, this changes the normalization of 
eq. {TJ: 



Ll ~ 30 I 1.4 



0.78 



0.6 7i, r 0.78 



L l» b M 



(20) 



Consequently, the normalization of eqs. (18191 changes as well: 



xs 



jet 



30 



1.4 

To 



-P 



dN 



BH 



(21) 



jet 



The numerical evaluation for j3 = 1.6 in the case of HMXBs 
gives only a slight change of 1.1, while [3 — 1.4 in the case 
of LMXBs gives 1.4. Note, however, that both £ and the 
mass are going to be different for neutron stars. Taking these 
into account, the absolute normalization the kinetic luminos- 
ity function is increased by factors of 2.4 and 3.0 for neutron 
star HMXB and LMXB respectively, compared to black hole 
H MXBs and LM XBs. 

iFended d2005 ) suggest that high-field neutron star systems 
do not produce jets, as they show no radio emission. Since a 
large fr action of the neu tron star HMXBs are, in fact, X-ray 
pulsars llLiu et al.ll200 0) with associated high fields, the nor- 
malization of the kinetic luminosity function for HMXB neu- 
tron stars is likely smaller by at least a factor of 2 compared to 
the above estimate. Furthermore, all HMXBs are believed to 
have rather strong fields (B > 10 12 G), which lFenderl J2005I) 
conjectured as being unable to produce jets at all. We there- 
fore conservatively assume that the HMXB values below are 
upper limits. 

Taking into account the lower X-ray luminosity that cor- 
responds to 77i cr it = 0.01 in neutron stars, the esti- 
mates for the mean kinetic power from neutron star jets is 
(W)hmxb < 2.6 x 10 37 ergs s^W^.s and (W) LM xb ~ 
2.5 x 10 37 ergs s _1 W37.8. We can then estimate the to- 
tal kinetic power from neutron stars into the Galaxy to 
be Wlmxb < 3.4 x 10 38 ergs s^Ws? 8 for HMXBs and 
Wlmxb ~ 3.3 x 10 38 ergs s^Wst.s for LMXBs. This brings 
the total kinetic power output from XRBs (black holes and 
neutron stars) to Wxrb ~ 9 x 10 38 ergs s -1 . Finally, the total 
energy released by neutron star LMXB jets over the lifetime 
of the Galaxy is -Elmxb ~ 2.2 x 10 56 ergsWsr.s^io- 



4. A COROLLARY: THE RADIO LUMINOSITY FUNCTION AND 
RADIO LOGN-LOGS OF XRBS 

Given the X-ray luminosity function for Galactic XRBs in 
eqs. and (|6) and the FP relation from eq. Q, we can easily 
write down the predicted Galactic radio luminosity function 
for this population of binaries: 

dN dN dl x 

— = -j-j--^- (22) 



where, following the nomenclature from above, l r = 
Lr/LEdd is the radio luminosity in units of the Eddington lu- 
minosity of a one solar mass object. The top axis in Fig. ^ 
shows the predicted cumulative radio luminosity function. 

For Galaxies with stellar populations of a similar age to the 
Milky way, the LMXB distribution should just be proportional 
to the mass M, while the HMXB distribution should be pro- 
portional to the star formation rate M: 

rfiV HMXB =10 _ 8 . 7/ _ 2M o.78 Col (23) 



dl T 

dNi^MXB 
dlr 



= io-°-'i 



5.2t-5/3 jiy-0.52 



Mire 



0.1 



(24) 



for 10" 9 < l x < 3.2 x 10~ 7 m° r f t . Again, while the total 
emitted radio power in the HMXB distribution seems to di- 
verge logarithmically if the lower limit i r ,min g° es to zero, 
the luminosity function is, of course, limited by the number 
of sources, which is not infinite. Furthermore, the radio spec- 
trum will eventually become optically thin at low luminosities 
and the luminosity function will become flatter than the ex- 
pression in eq. J22l > and converge. Note that in X-ray binaries 
this transition happens at very low luminosities: For typical 
low/hard state sources, the break from optically thick to thin 
occurs in the infra-red, and the break frequency v\> roughly 

l20?)l . 



follows fb OC Lr^ 17 



(Heinz & Sunyaev 2003). For the break 
to move below radio frequencies, we would have to consider 
radio fluxes 8 orders of magnitude below what is typically ob- 
served. 

Finally, we can estimate the flux distribution from the radio 
luminosity function (log N - log S). The flat spectrum XRB 
radio flux distribution for a galaxy at distance D is simply 
d log N/d log S — 4irD 2 dlogN/d\ogl I , since all sources 
are essentially at the same distance. Unfortunately, it is clear 
from the range in luminosities implied by Fig.^that sources 
in Galaxies further than about a Mpc will not be observable 
any time soon unless the population is large enough to contain 
significantly beamed sources (note that the verdict on how rel- 
ativistic XRB jets actual ly are is s t ill out - see, for ex ample, 
the discussio n inlGallo et alJj2003l):lHeinz & MerloniH2004l, : 
Kaiser et al] fc004 : lNaravan & McClintock H2005l) '). 

To derive the predicted radio flux distribution in the Milky 
Way, we have to convolve the luminosity function with the 
space distribution of Binaries inside the Milky Way. Follow- 
ing ( Grimr rTet alJ2002l) and taking r and z to be the radial and 
vertical distance to the Galactic center, the disk populations is 
describe by the function 



Tldisk^, z) 



(_Hn_JL_ ill) 

e r " d r * 8irr m rdr z 



(25) 



with r m w 4 kpc, r d w 3.5 kpc, r z w 410 pc for LMXBs 
and r m w 6.5 kpc, rd w 3.5 kpc, r z « 150 pc for HMXBs, 
and K2 is BesselK. The LMXB distribution also requires a 
spherical halo component, which we will assume to follow 

nh ^ R) = i^wm [kj (26) 

with b w 7.7 and R c w 2.8 kpc. Finally, about 30% of the 
disk sources reside in the Galactic center population, at about 
8 kpc distance, for which we use the distribution 



^bulgc 

(R) 



27rr[0.4] R^Rl 2 



(27) 
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FIG. 2. — Estimate of the radio flux distribution of Milky Way X-ray bi- 
naries observed at solar radius. Grey: HMXBs, black: LMXBs. Thick solid 
lines: Z m ; n = 10~ 4 ,'m cr ; t = 0.01, thin dashed lines: i m ; n = 0,rh cr ; t = 
0.01, thin solid lines l min = 10 -4 , merit = 0.001. 



with R = 1.0 kpc and R t — 1.9 kpc. 

The fraction of LMXBs in the halo is approximately 25% 
(including Globular Cluster sources). We will take the func- 
tional form of the luminosity function to be the same for the 
halo, bulge, and disk components. We then integrate over the 
luminosity function plotted in Fig.^and arrive at the flux dis- 
tribution shown in Fig. [2] for different values of upper and 
lower cutoff Z m ; n and m cr j t . The thick lines indicate the fidu- 
cial luminosity interval of 10~ 4 < l x < 10 _1 . Note that 
these are time averaged curves. Temporal variability neces- 
sarily introduces large uncertainty at the high flux end, where 
few sources contribute at a given time. 

Since we know at least one HMXB source which emits reg- 
ularly at 10 mJy levels (Cygnus X-l) and a spectrum of other 
transient sources that regularly pass above the lOmJy line, the 
HMXB curve derived for the fiducial parameters can safely 
be regarded as a lower limit. 

Clearly, the question of how far the X-ray luminosity func- 
tion extends to lower luminosities and where the radio-X-ray 
relation breaks down has large implications for the number 
count of sources at lower fluxes. In other words, extending 



the sensitivity of XRB monitoring campaigns to lower fluxes 
will reveal rather quickly below which radio luminosity the 
predicted radio luminosity function breaks down. 

Finally, it should be noted that this is essentially the prob- 
ability distribution of finding a given binary at a given dis- 
tance from earth, convolved with the predicted radio luminos- 
ity function. Since the actual number of XRBs is small, the 
error bars on this curve especially for large fluxes are proba- 
bly large. However, a realistic assessment of the uncertainties 
in this distribution is beyond the scope of this paper. 

5. CONCLUSIONS 

Starting from the observed radio-X-ray-mass relation for 
accreting black holes and the observed X-ray luminosity func- 
tion for XRBs, we derived the predicted kinetic luminosity 
function of compact jets from Galactic black holes in the 
low/hard state. The integration over the kinetic luminosity 
function of compact jets yields estimates of the average ki- 
netic power output of (Wjet) ~ 2 x 10 37 ergss _1 (dominated 
by LMXBs) and total integrated kinetic energy input over the 
history of the Galaxy of £tot ~ 4.5 x 10 56 ergs. We argue 
that transient jets should carry a comparable amount of power. 
We also derived the predicted radio luminosity function and 
the radio flux distribution for XRB jets. These estimates can 
be used in future modeling of XRB jet parameters and pro- 
vide a base line for estimating the impact of XRB jets in the 
interstellar medium. 
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